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ABSTRACT 

Theories of thick disk formation can be differentiated by measurements of stellar elemental 
abundances. We have undertaken a study of metal-poor stars selected from the RAVE spec- 
troscopic survey of bright stars to establish whether or not there is a significant population of 
metal-poor thick-disk stars ([Fe/H] — 1.0) and to measure their elemental abundances. In this 
paper, we present abundances of four a-elements (Mg, Si, Ca, Ti) and iron for a subsample of 212 
RGB and 31 RC/HB stars from this study. We find that the [a/Fe] ratios are enhanced implying 
that enrichment proceeded by purely core-collapse supernovae. This requires that star formation 
in each star forming region had a short duration. The relative lack of scatter in the [a/Fe] ratios 
implies good mixing in the ISM prior to star formation. In addition, the ratios resemble that of 
the halo, indicating that the halo and thick disk share a similar massive star IMF. We conclude 
that the a-enhancemcnt of the metal-poor thick disk implies that direct accretion of stars from 
dwarf galaxies similar to surviving dwarf galaxies today did not play a major role in the formation 
of the thick disk. 
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Introduction 



Thick disks are a common and significant stellar component in most disk galaxies, and their formation 
is an integral part of disk-galaxy formation. Most of the sta rs in the thick disk of the Milky Way Galaxy 
are old, ^,10 Gyr |Gilmore fe Wvse 1985 : Reddv et al. 20061) . and so they can serve as fossil records of the 
formation processes in early Galactic evolution. An important method for unlocking this information is the 
analysis of the elemental abundance patterns in thick-disk stars. Of particular importance are the ratios of 
the a-elemental (eg. Mg, Si, Ca, Ti) abund ances to iro n, which provide information about the past star 
formation and IMF of a stellar population (cf. IWvsdl2010l ). This approach is complemen tary to comparisons 
of the age distribution of thick-disk stars with theoretical expectations (e.g. I Wvsd 120091 ) . 



Model s of the formation of the t hick disk, including scenarios ranging from migrati on of stars from the 
inner disk ( Schonrich fe Binnev 20091) to heating of the thin disk due to mergers (e.g. Villalobos fe Helmi 
2008), make specific predictions about the chemical abundance properties of the metal- weak (and oldest) 
stellar population in the thick disk. In hierarchical clustering, heating of pre-existing thin stellar disks by 
merging continues until late times. In this case, an intermediate-age thick disk would emerge from heating 
of the thin disk, in conflict with the old ages found. A prop osed solution is tha t the old thick disk is a result 
of direct accretion of old stars from a few satellite galaxies (jAbadi et al.ll2003l ). The metal-poor stars in the 
thick disk would then show the same chemical enrichment as the parent satellite. Depending on the mass, 
orbit, and density profile of each accreted satellite, the a-abundance patterns of the thick disk would also 
vary with Galactic radius. The identification and chemical analysis of metal-poor thick-disk stars beyond 
the solar neighborhood is therefore very important. 



The existence of a significant metal- weak ([Fe/Hl < 



poor thick disk was identifi ed in both the field (jNorris 



1995 



et al 



1.0) thick disk remains controversial. The metal- 
Morrison et al. 1990t Wvse fc Gilmorel 



1985 



Chiba fc Beersl 120001 ) . and in globular clusters (jDinescu et al.l Il999i ). The apparent lack of blue 



horiz ontal-branch and RR Lyrae stars in the thick disk, however, argues against an old me tal-poor compo- 



nent ( Kinman et al 



Bensbv et all 12003 : 



2009). The small numb e r of metal-poor thick-disk stars studied to da te (jFulbri ght 2002; 



Brewer fc Carnevl 120061; iReddv et all bood Eeddv fc Lambertl boosl ) have high [a /Ye 



ratios, implying they formed in a short-duration star formation event, and old ages when estimates are 
available. These samples, however, have been limited to very bright stars within the solar neighborhood. 



The unprecedented size of the Radial Velocity Experiment spectroscopic survey (RAVE. ISteinmetz et al 

2006) and its selection without either kinematic or metallicity criteria provide a unique opportunity to study 
a statistically significant sample of bright metal-poor thick-disk stars, as described below. Our sample both 
probes a much larger volume and extends to lower metallicities than any previous thick-disk sample with 
elemental abundances. We here present elemental abundance results for a subsample of (predominantly) red 
giant branch (RGB) and red clump/horizontal branch (RC/HB) stars. 
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2. Observations 

2.1. Candidate Selection 

Candidate metal-poor thick-disk stars were selected from RAVE, a magnitude-limited survey that uses 
the 6dF spectrograph on the UK Schmidt telescope to obtain high S/N, TZ ~ 7500 spectra (A8410 A-A8795 A) 
of stars in the southern sky with / < 13. The RAVE pipeline provides estimates of radial velocities and 
stellar parameters (effec tive temperature, T c s\ gravity, log g; metallicity, [M/H]) through fits to a grid 
of synthetic spectra (see Zwitter et al. 20081 ). The stars are bright enough that they have proper motion 



measurements in the literature (included in the RAVE database). We selected candidate metal-poor disk 
stars with [M/H] < —0.7 and estimates of 3-D space motions, based on parameter values in the database, 
consistent with disk kinematics. 



2.2. High Resolution Echelle Observations 

High resolution spectroscopy provides robust abundances for most elements. We obtained data be- 
tween May 2007 and February 2009 using the following spectrographs/telescopes: MIKE/Magellan-Clay, 
FEROS/MPG 2.2-m, UCLES/AAT, and ARCES/APO-3.5. All instruments have a resolving power be- 
tween 35,000-45,000 and a spectral coverage of A3500 A - A9500 A, except for UCLES which covers A4460A- 
A7270A. The raw FEROS spectra were reduced using the Data Reduction System within ESO-MIDAS, while 
all other spectra were reduced using the echelle package in IRAF0. The final spectra yielded a S/N > 100 
per pixel at A5000 A - A6000 A and a minimum S/N ~ 40 around A4000 A, sufficient for abundance analysis. 

A total of ~ 500 spectra for candidate metal-poor thick-disk stars were obtained. We removed potentially 
problematic stars (e.g. hot stars, fast rotators), in addition to 4 radial- velocity outliers, probable binary stars, 
after comparison with RAVE (AT^ > 20 km s _1 ). Ten of the remaining stars have repeat echelle observations. 
We here report on the subsample of evolved stars (log g < 3.3), consisting of 212 RGB stars and 31 RC/HB 
stars, for which a uniform analysis is implemented. 



3. Abundance Analysis 



Our initial an alysis followed th e methodology of Fulbright ( 2000. hereafter FOOh . We combined the line 
lists from F00 and I Johnson! (|2002h . thus extending t he analysis to ext remely metal-poor stars. Equivalent 
widths (EWs) were measured using the ARES code (jSousa et al.1 120071 ). which fits a Gaussian to each line. 
Strong lines with EW > 110 mA were removed to minimize poor fits. Comparisons with hand-measured 
EWs for 12 metal-poor stars resulted in a mean difference EWhand — EWares ~ —0.1 ± 3 mA, sufficiently 
small that line-measurement biases are unimportant. 

The MOOG analysis program (|Sneder]ll973l) was utilized in an iterative procedure to compute elemental 
abundances using 1-D, LTE, plane-parallel Kurucz model atmospheres H. The stellar temperature was set by 
the excitation temperature method based on Fe I lines. The log g value was set by minimizing the difference 
between the calculated abundance of iron from the Fe I and Fe n lines. A microturbulent velocity was selected 



1 Distributed by NOAO, operated by AURA under cooperative agreement with the NSF. 
2 http:/ /kurucz. harvard.edu/. 
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to minimize the slope of the relationship between the iron abundance derived from Fe I lines and the value 
of the reduced width of the line. The metallicity value of the stellar atmosphere was chosen to match the 
iron abundance in the analysis. The repeat observations gave internal errors of 55 K, 0.1 dex, and 0.06 dex 
in T e g, log g, and [Fe/H], respectively, consistent with those of F00. 

External errors are critical in distance determination. We tested the accuracy of our parameters using 
echelle data (obtained for another project) of several globular cluster RGB stars (log g < 2) for which the 
distances and metallicities are reported in the literature. We also reanalyzed a sample (herea fter FOO/Hip) 
of 6 g iant stars from F00 which had Hipparcos parallaxes with acceptable errors (c p /p < 0.2; Ivan Leeuwenl 
20071 ). We found that with our initial parameter values many globular cluster giants lay above the RGB-tip, 
contrary to their location in the CMD, implying that our initial log g estimate is too low. Further, the 
interplay between log g and T e g in our initial analysis implies that our T c ff estimate would also be suspect. 
We therefore computed independent estimates of both gravity and temperature to quantify these effects and 
hence correct them. 

The independent estimate of surface gravity for each globular cluster and FOO/Hip star was obstaincd 
using the equation <7 DO i = AirGMa T* s /L. We adopt our initial estimate of T e g, a typical RGB star mass of 
0.8 Mq , and estimated the luminosity f r om th e de-reddened 2MASS K§ with bolometric corrections derived 
from Gonzalez Hernandez fe Bonifacio ( 2009^ and the published distance. We also derived an independent 
photometric temperature, TnViot,, for ea ch star using the 2MASS color-temperature transformations from 
Gonzalez Hernandez fe Bonifacio! (|2009l ). We assumed that the scale provided by T p hot does not show spu- 
rious trends such as apparently introduced by our initial spectroscopic analysis, and thus can be used to 
correct our log g estimates (as achieved below). 

Comparisons between our spectroscopic estimates and these derived values revealed offsets in both 
temperature and gravity. Further, the difference between our spectroscopic T e s estimates and T p hot is 
correlated with iron abun dan ce such that our estimate is 300 — 400 K cooler for the lowest metallicity stars 
(note both Johnson 20021 and Aoki et al. 2005 found similar offsets, albeit using different photometric colors 
in the derivation of T p hot)- We therefore corrected T e ff according to this correlation and repeated the analysis 
to obtain a new (ionization-balanced) gravity estimate, log g p hot ■ Stars with log <? p hot > 1 showed no mean 
offset with log g^oi and we therefore adopted log <? p hot as our final gravity. The difference, log <? p hot — log gbol 
is however correlated with log <? p hot for log <? p hot < 1- For these stars, we adopt a new log g estimate from 
this correlation along with the corrected temperature above to get a final estimate of the iron abundance. 
The final temperatures and gravities showed scatter with respect to Tphot and log gboi of 140 K and 0.2 dex, 
respectively. In the cases when the iron abundance from Fei and Fen do not agree, we chose Fen as our 
final estimate, since Fe II is both the dominan t species and much less sensitive to non-LTE effects than Fe I 
(jThevenin fe IdiartJ[l999[ Usplund et al.lll999h . The iron abundances derived for any given globular cluster 
showed a star-to-star scatter of ±0.1 dex, consistent with our earlier estimates of our internal errors. 



We used the above procedure to insure that the derived parameter values for our RAVE sample provide 
accurate distances. 
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4. Sample Analysis 
4.1. Distances 

Distances are critical and a full description of our technique is given in Ruchti et al. (2010, in prep). 
In short, di stances to our stars are estimated from the most probable Mx s magnitude in a grid of Padova 



isochrones ( Marigo et al. 2008t Girardi et al. 2002 ). The temperature, gravity, and metallicity, following the 



procedure in §3, and their associated errors (as described above) of 140 K, 0.2 dex, and 0.1 dex in T c g , log g, 
and [Fe/H], respectively, are used in the multi-parameter fit. For the RGB stars, our technique weights each 
point by the lu minosity function, derive d from the BaSTI luminosity function tracks (these isolate the RGB 
from the AGB; iPietrinferni et al.l 120041 ) . to emphasize evolutionary stages with longer lifetimes. Several of 
our stars have T e g and log g inconsistent with the RGB, and appear to be RC/HB stars. The position 
of a star on the horizontal branch depends on the mass- loss on the RGB, which is not well understood or 
modeled. The Padova isochrones, which we use, represent the zero-age HB as one point. Stars within 2<7T off 
of this point were assumed to have an absolute magnitude equal to that of that point. Note that Mk s ranges 
between only 0.1 — 0.2 magnitudes during the slow-evolution phase of core-helium burning. We assumed 
an old age for the RGB stars, performing a weighted average over isochrones of ages 10, 11, and 12 Gyr 
(reducing the age to 5 Gyr increased the distances by only 10%), while we assume an age of 12 Gyr for the 
RC/HB stars. 

We first applied this technique to the cluster and FOO/Hip stars. Our distances differed from literature 
values by only —2 ± 15% for the cluster stars and —4 ± 13% for the FOO/Hip stars. We therefore adopted 
a conservative estimate of 20% error on the distance, including both sc atter and offset. Distance estimates 
based on RAVE pipeline values of stellar parameters are now available! Breddels et al. 201ot Zwitter et al. 
2010(). Our distance estimates for the 172 stars with | log g e cheiie — log 5rave| < 0.5 are shorter than those 
of Zwitter et al. by 15 ± 24%. Note that our technique was optimized for metal-poor stars with echelle 



parameters, while IZwitter et al.l (|2010l ) optimized their method for all stars in the RAVE catalog, which 
have a high mean metallicity and typically younger ages. 

The average distance to our RAVE stars is ~ 2 kpc; all within ~ 7 kpc, except one at ~ 16 kpc. The 
majority of our stars have Galactic longitudes between t = 200° and 50° and Galactic latitudes \b\ > 25°, 
extending to an average vertical height of \z\ ^ 1 kpc. 



4.2. Population Assignments 

The 3-D, cylindrical space motion (Vn, Ve, Vz) of each star was computed by combining the de- 
rived distanc es and radial velocities w ith proper motions from the RAVE database. We adopted (10, 5.25, 
7.17) km s~ 1 ( Dehnen fc Binnev 1998h for the solar motion^ and Vlsr = 220 km s _1 . Each component of a 
star's space motion was then re-sampled 10,000 times, assuming a normal distribution centered on our esti- 
mate of the component velocity, with standard deviation equal to the propagated error in the velocity. We 
computed the probability that each re-sampled value of the space motion was drawn from a given kinematic 
population, based on the combined local characteristic Gaussian distributions assumed for each Galactic 



3 Our assignment method is insensitive to values within the range of recent estimates (e.g. Schonrich et al.ll201Clh . 
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population (see Table [T]), denned as: 



P(V n ,V e ,Vz)txexp --\ - 



r 2 



2a n 



(Ve - 220 - (V e )) 2 
2^1 



24 



(i) 



where, for each population, an, ce, and oz are the characteristic velocity dispersions and (Ve) is the 
mean rotational velocity (corrected for the Sun's motion). Note that Gaussian di stributions ar e a first-order 
approximation, and the distribution functions may be much more complex (cf., Binnev 2010[) . We did not 
include vertical gradients in kinematics or metallicity-kincmatical relations at this point. 



Table 1: Local Characteristic Velocity Distributions. 





on 




oz 


(Ve) 


Ref. 






(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 






Thin Disk 


39 


20 


20 


-15 


Soubiran et al. 


(2003) 


Thick Disk 


G3 


39 


39 


-51 


Soubiran et al. 


(2003) 


Halo 


141 


106 


94 


-220 


Chiba & Beers 


(2000) 



For each of the 10,000 samplings of the space motion of a given star, an assignment to the thick disk 
was made if that probability was 4 times that of either the thin disk or halo, with an analogous procedure 
for thin disk and halo. Intermediate populations, thin/ thick and thick/halo, were assigned for lower values 
of probability ratios. A star was assigned to the Galactic population with the highest number of occurrences 
from the 10,000 re-sampled points. 

A second population assignment was determined by comparing a star's position i n the Galaxy to the 
characteristic density distributions of each Galactic component (from ljuric et al.ll2008h . following the same 
Monte Carlo procedure as the space motion criterion. This criterion works well far from the Galactic plane, 
however discrimination is difficult close to the plane since a star must have \z\ = 2 kpc before the thick disk 
probability is 4 times that of the thin disk. The assignment is therefore restricted to choosing only halo over 
thick disk or thick over thin disk. 

The final population assignment was then the result from the space motion criterion, unless overridden 
by the boundary condition set from the star's position. Our results are consistent with each star's position in 
the reduced proper motion diagram. These final assignments were 73 thick disk, 1 thin disk, 22 thick/thin, 
31 thick/halo, and 116 halo. The Toomre diagram (Figure [T]) illustrates the correspondence between these 
population assignments and velocities. It is not surprising that many candidates were assigned to the halo 
since, with our final parameter values, our sample extends to large distances and high velocities. 



5. Alpha Element Abundance Results 

Elemental abundances were derived through MOOG after all stellar parameter corrections in §3. Figure[2] 
displays [a/Fe] as a function of [Fe/H] for several a-elements. Note that the range of [Fe/H] shows that our 
selection function from RAVE was efficient. The main bulk of thick-disk stars extend to [Fe/H] ~ —1.8, 
with 5 additional stars assigned to the thick disk below -2 dex. Further, these stars have [a/Fe] significantly 
above solar, with relatively low scatter, and blend smoothly into the halo stars. 



The kinematics of the metal-poor thick disk may differ from the canonical thick disk (jCarollo et al. 

2010|), with higher velocity dispersions and a slower rotational velocity. Adopting (Ve) = —100 km s^ 1 (see 
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Fig. 1. — The Toomre Diagram for our sample with ov e < 100 km s _1 . The black circle, green diamonds, 
red squares, orange plus signs, and blue triangles correspond to thin disk, thin/thick, thick disk, thick/halo, 
and halo stars, respectively. The dashed curves indicate const ant space motion , in steps of 100 km s -1 . 
The thick solid curve is an estimate of the local escape velocity ([Smith et al.ll2007j) . The total velocity error 
depends on distance and proper motion errors, with the typical la error < 20 km s . The star near the 
estimated escape velocity is the subject of a future paper. 



Gilmore et al.ll2002l ). would reassign several stars with thick/halo population assignments to the thick disk. 
There would be no change in our conclusions, since these stars have similar a-enhancement to the metal-poor 
thick disk and halo. 

A two- component h alo at [Fe/H] Si — 1, separated in [Mg/Fe] enhancement, has been reported by 
Nissen fc Schuster ( 2010l ). Indeed, several of our halo st ars in this metallicity range have values of [Mg/Fe] ^ 0.2, 
but other a-elements show little spread, in contrast to Nissen fc Schuster ( 2010h . We will return to this in 
another paper. 



6. Discussion 



We have computed the abundances of four major a-elements for RGB and RC/HB stars in our sample, 
with —0.5 ^ [Fe/H] }t — 2.8. About 40% of these stars have the highest probability of being (thick or thin) 
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disk stars and probe distance s much further than any earlier investigation. Previous high-resolution samples 
(e.g. Reddv fc Lambert]|2008f) did identify metal-poor thick disk stars with [Fe/H] < —1.5. We have increased 



this number by almost an order of magnitude, finding 24 thick disk stars with [Fe/H] < —1.5 that extend 
to metallicities below -2 dex. The [a/Fe] ratios for these stars are enhanced, indicating that the metal-poor 
thick disk was enriched by SNe II, before SNe la enrichment, implying a timescale for star formation shorter 
than ~ 1 Gyr. The relatively low scatter in [a/Fe] ratios also implies that the ISM was well mixed prior to 
star formation, and the metal-poor thick-disk stars were enriched by supernovae from an invariant massive 
star IMF. Further, the metal-poor thick disk and halo stars show similar a-enhancement, evidence that they 
were pre-enriched by the same massive star IMF. 

The a-enhancement in the metal-p oor thick disk contra sts with that seen in stars of similar metallicities 



in local dwarf galaxies (as reviewed bv lTolstov et al.1 120091 ) . All local dwarf galaxies began star formation 



earl y, 10 — 12 Gyr ago, and had extended star formation, often with dominant intermediate-age populations 



Tolstoy et al]|2009l). Dire ct acc retion of stars by assimilation of dwarf satellite galaxies into the thin and 



(e.g. 

thick disks in the lAbadi et al.l (|2003l ) models lasts until z ~ 0.7 (~ 6 Gyr ago). If the accreted dwarf galaxies 
formed stars until accretion or had extended star formation similar to surviving dwarfs, then many of the 
accreted stars will have formed from gas which had significant iron contribution from SNe la. Consequently. 



stars accreted from dwarfs into the thick disk after this time will have low [a/Fe] ratios ([Unavane et al 



1996), contrary to the enhancement we see in our sample. We conclude that direct accretion of stars from 
dwarf galaxies similar to surviving dwarf galaxies today did not play a major role in the formation of the 
thick disk. 

Another possibilit y is that the thick disk was formed from multiple gas-rich minor mergers, as simulated 
bv brook et all |2005l ). This model predicts a high star formation rate from the gas dissipation, consistent 



with our enhanced [a/Fe], but the model includes direct accretion of stars from the satellite galaxies. We 
would therefore still expect to find some low [a/Fe] stars in the thick disk from late merging, which conflicts 
with our low scatter in the abundance ratios of the metal-poor thick disk. Early heating of a thin stellar 
disk by mergers, however, is still viable. 

The amplitude of radial and vertical gradients in abundances and metallicity will provide diagnostics 
to further discriminate models of the formation of the thick disk. Our sample, which probes distances much 
further from the solar neighborhood, is the first for which this can be done, and will be investigated in our 
next paper (Ruchti et al. 2010, in prep). 
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Fig. 2. — Computed [a/Fe] ratios versus [Fe/H] for our sample of stars. While Fen is used to estimate 
[Fe/H], element ratios are computed using th e iron abundan c e of t he same ionization state as the a-element 
(e.g. [Si/Fe] = [Sii/Fei]) as is sugg ested bv iKraft fc Ivansl |2003h . Color and symbols are the same as in 
Fig. [T] The cause of the offset between [Tii/Fe] and [Tili/Fe] is unclear, but the thick disk and halo still 
show similar enhancement in each. Two stars of note are the halo star at [Fe/H] ~ —1.45 with very low 
a-enhancement and the thick disk star at [Fe/H] ~ — 1 with a very high [Si/Fe] ratio, each of which is the 
subject of future papers. 



